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K— 5 Abstract 

OS In this paper, we analyze the dynamics of an N particles system evolv- 

ing according the gradient flow of an energy functional. The particle sys- 
tem is a consistent approximation of the Lagrangian formulation of a one 

fjy parameter family of non-local drift-diffusion equations in one spatial di- 

mension. We shall prove the global in time existence of the trajectories 
of the particles (under a sufficient condition on the initial distribution) 
("^ and give two blow-up criteria. All these results are consequences of the 

competition between the discrete entropy and the discrete interaction en- 

Sergy. They are also consistent with the continuous setting, that in turn 
is a one dimension reformulation of the parabolic-elliptic Keller-Segel in 
high dimensions. 
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O 1 Introduction 

m 



X 



sive particles 



This paper aims to give an insight of the blow-up dynamics driven by the com- 
petitions between random and directed movements undergone by particle sys- 
tems. More specifically, we consider the following one parameter family of one 
dimensional non-local drift-diffusion equations describing self-attracting diffu- 



d t p = d xx p-d x (pd x S y (p)), xeM.OO. (1.1) 



Here p is the density of the particles and 5 7 (p) = K~ * p is the interaction 
potential, with the interaction potential kernel if~ given by 



K 7 (x) := 7 -1 |ie|- 7 , 76(0,1). (1.2) 



Non-local equations of type (1.1) in space dimension d > 1, have been in- 



troduced in several domains. Let us recall the parabolic-elliptic Keller-Segel 
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system [HO HOI E] , a first attempt to model the aggregation of Dyctiostelium 
discoideum, the spatially homogeneous and dissipative Boltzmann equations 
for (simplified) granular flows [21 O [13l [24] , the Smoluchowski-Poisson equation 
modeling gravitational collapse [H [5J [3] , and finally the McKean-Vlasov equa- 
tions in kinetic theory [24] . In any of such mathematical model, the spatial 
dimension and the interaction potential kernel are chosen accordingly to the 
physical or biological phenomenon under consideration. However, the interac- 
tion potential kernel is always assumed to be even, thus reflecting the Newton's 
third law. Furthermore, they all share the property to be endowed with a free en- 
ergy, given by the combination of an internal energy with an interaction energy. 
This common feature is due to the attempt to model random but self-attractive 
movements of particles, whose mass density is conserved along time. More in- 
terestingly, the free energy gives the variational formulation of the physical or 
biological phenomenon, since it is decreasing along the trajectory of the parti- 
cles and satisfies an "H-theorem" . As a consequences, the mathematical model 
can be seen (at least formally) as the gradient How of the energy in the space 
of probability measures endowed with the appropriate metric. The latter fact 
being rigorously proved nowadays for some of the previously cited models [U [7] . 

Despite all this common features, the analytical behavior of the solutions 
of those type of equations can be completely different according to the choice 
of the interaction potential kernel. More precisely, if the interactions between 
particles do not enter in competition with the diffusion process (which would be 



the case in ( 1.1 )-( 1.2 ) with a negative 7), then the combined effect of the internal 
energy and the interaction energy gives rise to global solutions converging toward 
the equilibrium [13l HU [24] . On the other hand, if the interaction potential 



kernel is in competition with the diffusion process, as in ( 1.1 1-( 1.2 ), then the 
solutions can blow-up in finite time whenever the initial density is "large" and/or 
"concentrated" enough. This is exactly the case of the parabolic-elliptic Keller- 
Segel system in dimension d > 2 ([8) [10]) and of the Smoluchowski-Poisson 
equation ([221 [23]) cited above. This kind of scenario is much less understood 
than the previous one, for which very accurate results exist (see the references 
in [H]). This is also the reason why we consider here the initial valued problem 



associated to (1.1)- (1.2) 



The family of equations ( 1.1 )-( 1.2 ) has no physical interpretation, to the 



best of our knowledge, but can be interpreted as the projection of the parabolic- 
elliptic Keller-Segel system from the high dimension d > 3 to the one spatial 
dimension. Indeed, the associated free energy is 

£y[p](t) = J^p(x,t)logp(x,t)dx- * J p(x, t) (A' 7 * p(t))(x) dx , (1.3) 

where the first term is the internal energy U (p) , also called entropy hereafter, and 
the second term is the interaction energy W 7 (p). Under regularity conditions it 
satisfies 



^£>](i)=- / p(x,t) \d x {log p{x,t)-(K 7 *p(t))(x))\ 2 dr. 
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Therefore, thanks to the decreasing behavior of the energy functional and ex- 
ploiting efficiently the competing convexities of the kernel K 7 and of the internal 
energy density p log p, it is proved that the initial valued problem associated to 



( 1.1 H 1-2 1 shows at least two blow-up criteria (see Proposition |2.2[), exactly as 



the high-dimensional parabolic-elliptic Keller-Segel system [H | 

However, the goal of this paper is not to analyzed the continuous model 
( 1.1 M 1.2 ), but a consistent approximation of its Lagrangian formulation. For 



that, taking advantage of the one dimensional setting, of the non-negativity of 
the solution p corresponding to a non-negative initial density p and of the 
conservation of the initial mass, i.e. 



p(x,i)dx= / po{x) dx =: M . 



t>0, 



(0,M) 



the family of equations ( 1.1 ) can be reformulated in term of the pseudo inverse 

X(;t) 



of the distribution function associated to p. For the 



pseudo-inverse function X, the energy functional (1.3 1 rewrites as 



gjx](t) 



M 



1 



log (d m X(m,t))dm-- 
o * Jo 



M r M 



K y (X(m, t)-X(m', t))dmdm! , 
(1.4) 

while the family of equations ( 1.1 )-( 1.2 ) rewrites as the integro-differential equa- 



tions (3.2). Moreover, the latter result to be the gradient flo w of the energy 
functional (1.4 1 in the Hilbert space L 2 (Q,M) (see Proposition 3.1). 



The particles approximation is then constructed from the mentioned gradi- 
ent flow interpretation, thereby preserving all the properties of the continuous 
problem. It results in the dynamical system (4.5) describing the trajectories 
X(t) = (Xi(t))i = x t ...,N of N > 3 particles distributed increasingly, i.e. lying 
in the set n = {X s M. : Xi < Xi+i}, and carrying the same fraction 
M(N+ 1) _1 of the conserved total mass M. For this dynamical system we give 



a global existence result (Theorem 4.1) based on a sufficient condition on the 



initial distribution of the particles that prevents collision. This sufficient condi- 
tion is far from being optimal, but it is in some sense better than the sufficient 
condition for global existence of weak solutions of the parabolic-elliptic Keller- 
Segel system in high dimensions. Moreover, we obtain two blow-up criteria for 



the particles system (Propositions |5 . 2| and 5.4) consistent with the blow-up cri- 
teria for the continuos equation (Proposition 2.2) and incompatible with the 



global existence result. More specifically, the particles whose initial distribution 
satisfies one of the blow-up criteria and the particles whose initial distribution 
satisfies the global existence condition, lye in the subsets of II complementary 



w.r.t. the curve of the critical points of the discrete energy (4.3). Along this 



curve the time derivative of 1^(^)1 is zero. However, the curve is not the un- 
stable manifold separating the two basins of attraction (global existence and 
blow-up), as it is put in evidence numerically in the case of a three particles 
system with zero center of mass. In that case, all the mentioned criteria have 
been plotted in Figure [T] and the unstable manifold computed numerically. The 
derivation of its equation would provide a single criterion to distinguish between 
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global existence and blow-up. It would be also the proof that the dichotomy 
of the limit case 7 = with logarithmic interaction kernel (see Section [6| still 
holds true in the case 7S (0, 1). The derivation of this single criterion is under 
investigation in a forthcoming paper. 

We conclude this introduction, observing that the Jensen inequality will 
play the main role in almost all the proofs. It is in fact quite the only tool used 
here. This reflects the fact that the competing convexities of the interaction 
potential kernel and of the internal energy density are the determinant keys 
of the dynamics of the particle system as well as of the continuous family of 
equations ( 1.1 1. 

This paper is organized as follows. In Section[2]we summarize the analytical 
properties of problem ( 1.1 )-( 1.2 ). In Section [3] we reformulate the problem in 
term of the pseudo-inverse of the distribution function associated to p and we 
give its gradient flow interpretation. In Section [4] we introduce the N particle 
scheme and we give a sufficient condition for the global existence of the tra- 
jectories. Section [5] is devoted to the analysis of the blow-up of the particles 
trajectories. Finally, in Section [6] we briefly recall the limit case 7 = and 
we give the corresponding global existence result for the trajectories of N > 3 
particles. 



2 Preliminaries 

In this section we briefly recall some facts and results about the existence and 
the blow-up of the solutions of the initial value problem for the family of equa- 



tions ( 1.1 )-( 1.2 ). These results are inspired from those obtained in the analysis 
of the classical parabolic-elliptic Keller-Segel system. They are given here for 
completeness, but also to convince the reader that the re- writing of the family 



of equations (1.1) in term of the pseudo inverse of the distribution function of 
p is the natural procedure to obtain explicit computations by its gradient flow 
interpretation, (see Secti on |3j ) . 



First of all, equation ( 1.1 ) is invariant under the space-time scaling p\(x, t) = 
A 7_1 p(^, w), which in turn conserve the L^fR) norm. This invariance prop- 
erty suggest that the Lebesgue space L 1 ^ (M) could be a critical space for the 
existence of global solutions of ( |l.l[ )-( |l~2"| ), as it is the case with hi (M. d ) for the 
high dimensional Keller-Segel system [10, US]. However, due to the singularity 
of the drift term d x S^(p), inherited by the strong singularity of the interac- 
tion potential kernel K~ ( in one dimension, this seems not to be the case for 



equation (1.1), (see also the end of Section [4]). More specifically, the evolution 



equation of the L p norm of p is given by 



i I P P dx = - A ^±\\d x pi\\ 2 L2(VL) +2{p-l) f p P W xP HK 7 *d x p)dx. (2.1) 
at Jr. P Jr 



From the weak Young inequality, the convolution with K 7 is a bounded map 
from L m (R) to L r (M), with 1 < m < r < 00 and £ = ±- + (7 - 1). Therefore, 
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for 1 = \ + § 



i + ( 7 -l),p6(l,2] andi = i 



(2-P) 



we have 



J p% d x p§ (K 7 * d x p) dx < C\\p*\\ L * 
9 

\d x p'< 



<-c||^IU, ( 
p 



2 

I iff 



L 2 -p ( 



li 2 (R) l|c?a;Plli"»(]R) 

I) p rll il-1l 



JIM 



2 

li=f 



(2.2) 



since, by the previous choice of indeces, ^ 



\2.l\ , it follows 



l 

1-7' 



Thus, plugging (2.2 1 into 



,/'(..-./),/,•< l^^||^l(t)||| 2(R) 



c\\ P {t)\\ 



L 1 — < ( 



1 



(2.3) 



dt .1 p 
However, we can not deduce any control of J R p 1 ^ dx from the previous differ 



ential inequality, since we are not allowed to chose p = in (2.3 ). Indeed, on 
the one hand 1 < p < 2, and on the other hand, by the previous numerology it 
holds: ~ > 0, implying q > 2, which in turn implies p < j^— . Consequently, it 
seems not possible to obtain a global existence result from a smallness condition 
on the L T ^> norm of the initial density pq. 



The previous computation implies that a weak formulation of problem ( 1.1 1 
should be considered. Owing to the property K 1 {x) — K 7 (—x), it is natural to 
consider 



dt 



p(x,t) 4>(x) dx = / p(x, t)(f>"(x) dx 



— y) 

((f)' (x) - </) (y))p(x,t) TZ — —^p(y,t)dxdy, 



(2.4) 



for any twice differentiable test function <p with bounded second derivative. 

Whenever p(x, 0) = po(x) in the sense of distribution, this formulation guar- 
antees that the total mass M is conserved. The existence of such a weak solution 
is an open question. A possible way would be to regularize the problem, as done 
for instance in j^T]. An alternative proof of existence could be obtained using 
the Jordan-Kinderlehrer-Otto scheme [HUH], as in [7]. However, we don't go 
through this matter since this is not the goal of the present paper. 



The weak formulation (2.4) allows also to obtain informations on "moments" 



of p. For instance, taking <f)(x 

of any weak solutions is conserved. Furthermore, for cj)(x 
evolution equation of the second moment lit) :— J R \x\ 2 p(x,t) dx of p, i.e. 



x in (|2.4|), it follows that the center of mass 

x 2 , it follows the 



I'(t) = 2M- p(x,t) 



p(y, t) dx dy . 



(2.5) 



nomenon [U E] (see also Section 



When 7 = 0, equation (2.5| gives the well known mass threshold phc- 

i 6|. In contrast with this simple dichotomy, in 
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the case 7 G (0, 1) considered here, at least two blow-up criteria for the density 
p can be deduced from (2.5). Indeed, (|2.5|) reads as 

(2.6) 



/' (t) = 2M - 7 W 7 (p) = 2M + 2 1 £ 1 [p]-2 1 U(p). 



It is then sufficient to obtain appropriate lower bounds for U(p) or W 7 (p) w.r.t. 
/ in order to get a differential inequality from (2.6) assuring blow-up. 

Lemma 2.1. Let 7 G (0, 1) and p G L^(K; (1 + |x| 2 ) dx). TVie following lower 
bounds for the interaction energy W 7 (p) and the entropy U{p) hold true 



7 W 7 (p) > 2-~</ 2 M 2+ i r~</ 2 , 



M 

U{p) > -y log/ 



M . /M 3 \ 

T log U^J 



(2.7) 
(2.8) 



Proposition 2.2 (Continuous blow-up criteria). Let 7 6 (0, 1) and assume that 
the initial density p G ^^(IR; (1 + \x\ 2 ) dx) has finite energy and satisfies either 



\x\ 2 po{x) dx < 



M 



or 



\x\ 2 po(x) dx < 



■ exp 



(2.9) 



(2.10) 



2ne^ +1 V M 

Then, if the corresponding weak solution has finite energy, the solution blow-up 
in finite time. 



A proof of (2.7) and (2.8) can be found in [4] and [10] respectively. However, 
we shall show in Section [3j an alternative proof that make use of the change of 
variable introduced there. Let us just observe here that the l.h.s. of (2.7) and 



The proof of Proposition 2.2 is also standard and not 



( |2^8| ) might be +00. 
given here. We only point out that both criteria (2.9) and (2.10) are invariant 
under the space dilation f\(x) = A 7_1 /(|) induced by the space-time scaling 
previously introduced. The interested reader can found the blow-up criteria for 



the high-dimensional parabolic-elliptic Keller-Segel system equivalent to (2.9) 



and (2.10) in [10 



3 The continuous equation as a gradient flow 

We shall hereafter take advantage of the one dimensional setting of the problem 
to perform an explicit "change of variable" allowing us to rewrite the family 



of equations (1.1) in a sort of Lagrangian formulation. For this purpose, we 



consider the distribution function associated to p, 



F(x,t) 



P(y, t) dy , 



x G 
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and its pseudo inverse X(-,t) : (0, M) —> R, 

X(m, t) := inf {x E R \ F(x, t) > m } . 

Since /?(•, t) E Li_(K.), F(-, t) is an absolutely continuous non decreasing function 
and therefore a.e. derivable over R, while X is a right continuous non decreasing 
function and a.e. derivable over (0, M) with d m X > (being X the gradient of 
a convex function, see [E]). Moreover, the following identity holds true for any 

m E (0, M) 

F{X(m,t),t) =m, (3.1) 

(since the measure pdx do not charge points, otherwise F(X(m,t),t) > m). 

Next, let us proceed with formal computations, without take into account 
the singularity of the drift term in (1.1 1. Integrating the latter over (— oo,x), 
one obtains the equation for F as 

d t F = d xx F-d x Fd x S 1 {d x F). 



Using the differential relations obtained from the identity (3.1 1, the previous 
equation reads as the following family of integro-differential equations 



d t X(m,t) = d„ 



1 



M 



X(m,t) - X(m',t) 
d m X(m,t)J ' J Q \X(m,t) - X(m> ,ip+ 2 



dm. (3.2) 



Obviously, the singularity in ( |1.1| has been translated into (3.2), giving rise to 
an integral term a priori not finite, and one should consider the following weak 
formulation 



d 
dt 



M 
1 

2 



(j){X{m, t)) dm 



M 



<f>"{X(m, t)) dm 



M i-M 



(4>'{X{m, t)) - 4>'{X(m', t))) K'(X(m, t) - X(m', t)) dm dm', 



(3.3) 



obtained taking into account the boundary conditions : 



<9 m X(0, t) = d m X(M, t) = +oo , X(0, t) = -oo and X(M, t) = +oo . 

(3-4) 

The previous results can be made rigorous, recalling that the map X trans- 
ports C\o t M]i the Lebesgue measure in the interval [0, M], over the measure pdx, 
i.e. X#C< 0> M] has distribution function F. With this transport point of view 
in mind, it is easily seen that (3.3) is the translation of ( |2.4[ ) in the m variable, 
that 



M 



X(m,t)dm= / xp(x,t)dx, 



and X(-,t) E L 2 (0,M) has soon has p has finite second moment, since 



X (m, t) dm 



\x\ p(x, t) dx . 



(3.5) 



(3.6) 
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Moreover, making use of the area formula, that can be applied here thanks to 
the regularity properties of the map X (see pQ page 130), the entropy U{p) 
becomes 



U(X)=- log (d m X) dm, 
Jo 



(3.7) 



and the energy (1.3) translates into (1.4) 



Despite of the singularity of the integro-differential equation ( 3.2 ) , we deduce 
finally the gradient flow interpretation of (|3.2|. 



Proposition 3.1 (Gradient flow interpretation). Let 7 g (0,1). The integro- 



differential equation (3.2) is the gradient flow of the energy functional Q^[X] for 
the Hilbertian structure over £ 2 (0, M): 



d t X = -VliG^X}. 



(3.8) 



Proof. To prove (3.8), we compute formally the first variation of the functional 
C? 7 LY], omitting the t variable for simplicity, to obtain 



A I 



d m Y{m) dm 



M i-M 



o d m X(m) 
K'^iXim^t) - X(m',t))(Y{m) - Y{m'))dmdm! 

n y/ s ) Y{m)dm 
d m X(m)) 

\ \ ^(X(m,t) - X(m',t))Y(m)dmdm' = - < d t X,Y > 
Jo Jo 







□ 

Remark 3.2 In order to obtain a rigorous proof of the previous proposition one 
should use the Jordan-Kinderlehrer-Otto scheme [Tl I18j. 

Remark 3.3 (Scaling invariance) The space-time scaling leaving invariant the 
equation ( |1.1| translates into X\(m,t) = \X(j^, ^-), the distribution function 
associated to p\(x,t) = A 7_1 p(|, ^) being F\{x,t) = \~ / F(f,j S ). Further- 
more, it is easily proved that the scaling X — > X\ leaves invariant the integro- 
differential equation (|3.2|. 



We conclude this section observing that (3.3) with (f>(x) = x 2 becomes (2.6), 
and giving the proof of Lemma 2.1 with the help of the transport map X. 



Proof of Lemma 2.1 We shall consider densities p with zero center of mass, 
without loss of generality. Inequality (2.7) follows immediately by the Jensen 
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inequality applied to the convex function <f>(z) = z 7 / 2 , and (3.5), since 

-M i-M 



7W 7 (p) 



p{x) 



1 



\x - y\ 



■ p(y) dx dy 



\X{m) - X(m')\-~< dmdm' 



> M M~ 



M r M 



JO 



\X(m) - X{m')\ 2 dmdm! I = (2I)~' M 2+ i . 



Next, let us assume M = 1, again without loss of generality. In terms of the 



pseudo inverse X, inequality (2.8) reads as 

log {X' (m)) dm +^log( J \X(m)\ 2 dm^ >i]QgfJL 



(3.9) 



Let p(x) := x G the distribution function associated to p and G 1 the 

inverse of G. Using G as a change of variable, the identities G"(x) = —xG'(x) 
and J m x 2 G'(x) dx = 1, the Jensen inequality again and the Cauchy-Scwartz 
inequality, we have 



log (X'(m)) dm 



and (3.9 1 is proved. 



log ((X(G(x)))') G'{x) dx + / log(G'(x))G'(x)dx 



> - log 
= - log 



(X(G(x)))'G'(x)dx ) - -log(27re) 



X{m)G- l {m) dm)-- log(27re) 



> -_ 



l^e(J \X(m)\ 2 drn^j -hog(2ne) 



□ 



Remark 3.4 It is worth noticing that inequality (3.9) cannot be derived using 
a different probability measure (because of the identity G"(x) = —x G'(x)), but 
it is invariant w.r.t. any dilation p\(x) = xp(t) of p. Moreover, it is an identity 
iff X — G^ 1 , i.e. p — p\, with A = ii , since 



X^G^im) dm 



\X{m)\ 2 dm+ i - i / \X(m) - G' 1 (m)\ 2 dm . 



4 The N particle scheme 



We are now ready to construct a finite dimensional dynamical system approxi- 



mating the intcgro-differential equation (3.2) for the pseudo inverse X, adapted 



to the gradient flow structure (3.8). The main point is to discretize the free 
energy £ 7 [-X] by using standard quadrature approximations and then to write 
the finite-dimensional gradient flow of the approximated energy. Following this 
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procedure, we guarantee the preservation of some important properties at the 
discrete level, e.g. the dissipation of the energy, the homogeneities of the diffu- 
sion and interaction contributions, the conservation of the center of mass and 
the competing convexities. 

To begin with, we introduce the regular grid of points = i h, i = 1, . . . , N, 
over (0, M) , where iV is a fixed integer greater than two (the number of particles) 
and h = M(N+ ly 1 is the space step (chosen constant for simplicity). Then, 
we denote by Xj(t) the approximation of X(m,-,t) and X(t) = (Xi(t))i=i,,„ t N- 
We make also the convention Xq = — oo and Xm+i = +00. That is to say we 
put ghost particles at ±00. This convention is consistent with the boundary 



conditions (3.4) and makes the computations below simpler. 

Since we approximate a nondecreasing function X(-,t) and because the free 
energy becomes singular whenever particles collide, it make sense to look for 
increasing trajectories X(t), i.e. lying in the set II := {X 6 R N : X t < X i+ i}. 



Then, for the entropy U{X) in (3.7) we opt for the approximation 



J\o g (d m x(m,t))dm « - h log { Xl+l{t \ Xl[t) ) =■■ mm), 

(4.1) 



while for the interaction energy W(X) we opt for the approximation 

i-M pM > 2 N N 

- / K 7 (X(m,t)-X(m',t))dmdm' w — V V K^{Xi{t)-Xj (t)) :=W y [X](t) 
2 Jo Jo 2 . =ij . =1 

(4-2) 

Therefore, we have derived as a natural approximation of C/ 7 [X] , the following 
discrete energy functional 

G^[X] := U[X] - W^[X] , XeU. (4.3) 

Accordingly, the finite dimensional dynamical system is defined as the gradient 
flow of G 1 in the euclidian space K. 



X{t) = --VG 7 [X](t). (4.4) 



1, 
h 

The dynamical system ( |4.4[ ) shares the same time/space invariance with the 
continuous setting (see Remark |3.3[ ) . It is in fact invariant under the rescal- 
ing (h\,X\(t)) = (A 7 h, \X(^)). Rescaling the step h is necessary in order 
to take into account that the mass M scales into A 7 M under the dilation 
f\(x) = A 7 ~ 1 /(|). Moreover, system (4.4) preserves the center of mass, i.e. 



YliLx Xi(t) = Si=i-^Q(0)- Indeed, using (4.1), (4.2) and the previous conven- 



tion, the dynamical system (4.4) reads in details as follows 



l{ ' X i+1 (t) - X t (t)^ Xi{t) - AVi(<) £j \Xj(t)- Xi(t)[r+i ' 1 ' 
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for i = 1, . . . , N. Then, by summing the differential equations (4.5 ) over i, using 
the telescopic summation, the boundary conditions, and the symmetry of the 
interaction kernel, we obtain Xi(t) = 0. 

The local in time existence of solutions X(t) £ H of the system of ODEs 
(4.5) is a consequence of the Cauchy-Lipschitz Theorem. However, the maximal 
time of existence can be finite. Clearly, this happens when at least two particles 
collide, or equivalently as the trajectory X{t) reaches the boundary of the set II. 
By the means of the function 



N-l 

E 



7 ^ (X i+1 - Xi)i 



we shall establish hereafter a sufficient condition on the size of the initial datum 
X(0) that prevents collision of particles, and assure therefore the global existence 
of the trajectory X(t). 

Theorem 4.1 (Global existence). Let 7 € (0, 1) and N > 3. There exists 
a constant c(N) > such that the following condition on the initial datum 

x(0) <= n, 

7 tfh(X (0)) < (1 + (N - 3)c(7V))- 1 , (4.6) 
guarantees the global existence of the solution X(t). 

Proof. We claim that (f>-y(X(t)) is strictly decreasing in time if X{t) is a tra- 
jectory starting from X(0) satisfying (4.6). Since <t> 1 (X(t)) goes to +00 as the 
trajectory X(t) reaches the boundary of the set II, the claim gives us the proof. 

Using (4.5), we can decompose the evolution of 7 along the flow X(t) into 
the diffusion contribution and the interaction contribution as following 



d 
dt 



N-l 

<£ 

i=l 



N 

i=l 



X i+l {t) - Xi{t) 
Xi{t) 



(X i+1 {t) - Xi(t))^ 



_{X i+1 (t) - Xi(t))7+1 

= -I 1 (X(t)) + J 1 {X(i)), 



(4.7) 



where 
J 7 (X) 



Xi+i — Xi Xi — Xi- 



{X l+1 - X t )^ (X t - Xi-ip+i 



is positive, and 



N N 



i=i j=i 1 3 



sgaiXj - X t ) 



X^+i \{X l+1 -Xi)^ (X, 
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In turn, J y (X) can be decomposed into the two interaction contributions due 
to contiguous and non-contiguous particles respectively, as following 

N / - ° 

j,(X) = h>J2 ' ' 



A N sgn(X j - Xi) ( 1 1 



^ 2 E E 



^ ^ \x, - V(x J+1 - xtf+i (x, - 

li-»l>2 
=: J 7 X (X) + J^(X) . 

(4.8) 

Moreover, with the notations K; := JQ+i — X,, i = 0, . . . , N, and Y = (K;),, for 
simplicity, we have the identity 



3=1 x i } \i=o \ Y i+i 



iV-l / , \ JV-1 



E ( v-y+1 V 7+l J E V 7+1 ( V 7+l V 7+l ) + 

(4.9) 



= h- 2 ji(Y) + R(Y) , 



with the positive remainder R(Y) 

f-r A; y? n+i I Kit 1 5-' +1 i f-: ^ K 1 « I Kit 1 k ,+1 



\ V~ +i V~ 

1 1 / 1 1 \ ^ 1 1 ( 1 1 ' 

V V ± (JL L\ (_J_ jJ\ 

+ l^i l^i vy \ y-.i y / \ W+ 1 VT+ 1 J 



v 1 v 1 jl r 1 m 

y? \ ylTi y / \ w +1 w +1 



Therefore, the evolution equation (4.7) becomes 

J^(*(t)) = iT^rM*)) - 1] Iy{X{t)) + J? t (X(t)) - h? R{X{t)) , (4.10) 
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and it remains to estimate the remainder J?(X(t)) — h 2 R(X(t)). The difficulty 
here is to measure the relative differences \Xj — Xi\, where \j — i\ > 2, w.r.t. 
\Xi±i — Xi\. For that reason we shall distinguish the cases of three and more 
particles. 

The three particles case. When N = 3, the above reminder is negative since 



h~ 2 J 2 {Y) - R(Y) = 



1 1 



1 



1 



1 



1 



1 1 



rl^i 1 " 7 -^" 7 ) 



l-7\ 



Y?Y? +2 Y?Y? +2 Y 1 Y 2 \YJ +1 Y[ 



7+ 



< 



{Y 1 Y 2 )i\Y 1 Y 2 J Y 1 Y 2 \Y? +1 Y[ 



7+ 



r M 



1-7 



The theorem is then easily proved in that case. 

The high number of particles case. When N > 3, we shall take advantage of 
the positivity of R(Y) and estimate J 2 (Y) w.r.t. J^(Y). Indeed, using the 
increasing behavior of the trajectory, telescopic summations and the boundary 
conditions Yq = — oo and Y/v = +oo, from 14.8b it follows 



N i-2 



i=3 j=l 



1 



1 



Y; 



7+1 



Y 



7+1 



N-2 N 



N - ( 1 1 

i=5 j=i \- I j-i z j 



i=l j=i+2 i 

1 1 

y7+l V7+1 
I i 



1 



1 



7+1 



N-2 i 
i=l 3=1 



7+1 



V7+1 
r 7-l 



1 



<(7V_ 2 )^ 



8=3 

AT-2 



^(F) + l(jV-i + l) 



y7+l Fj 7 ^ 1 
1 1 



7+1 



7+1 



i=i 



7+1 



7+1 



< 2(N ~2) 2 h- 2 JhY) 



Y : 



7+1 



(4.11) 



Then, plugging (4.11) into the differential equation (4.10) and using the identity 
(4.9), we obtain the differential inequality 

~£y(*(t)) < [(1 + - 2) 2 ) 1 ^{X{t)) - 1] I 7 (X(t)) , 
and the theorem follows again. □ 
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We conclude this section with some remarks on the smallness condition (4.6 1 
for global existence. The function 7 has been chosen for its homogeneity prop- 
erty, that allows us to obtain the differential inequality on </> 7 itself. In addition, 
4> 1 is invariant under the rescaling (h\,X\) = (A 7 h, X X) of the discrete setting. 
More interestingly, </> 7 (X) controls the interaction potential W~(X), since the 
following inequalities (also invariant) hold true for all X e II 

N 

h<t> 7 (X) < W y [X] < h — (j> 7 (X). (4.12) 

Therefore, it is natural to obtain a global existence criterion from the analysis 
of the evolution of </> 7 (A"(i)). Indeed, the discrete blow-up criteria in Section [5] 
are based, roughly speaking, on the fact that if the interaction energy Wj is 
initially large, then it remains large and aggregation dominates diffusion. 



The left inequality in (4.12) is an immediate consequence of the definition 
of Wj. For the right one, we have 



N-l N N j-1 

•2 \ " \ " 1 1 1,2 



i=l j=i+l v J *' j=2 i=l v J 11 



1 



N-l ^ N 



i=l ^ j=2 

N-l 



i=l 



Moreover, following the previous discretization procedure, it is easily seen that 
7/i 7 7 (X) is an approximation of IHIIrnnay Taking this into account, the 



right inequality in (4.12) results a (non-optimal) discretisation of the following 
inequality in the continuous setting (derived from the Hardy-Littlewood-Sobolev 
inequality plus interpolation) 

P( x ) rz 1 „,\ y P(y) dxd V ^ C '7 m1_7 IIpII2^+i(r) • 



The non optimality is due to the crude inequality established in ( 4.13| ). Since the 



L' t+1 norm is controlled by the presumably critical norm L 1 -^ , i.e. \\p\\'Jt+i(u) < 



^ 7 llp|l ir J— the smallness condition (4.6) on <f> 1 (X(Q)) is less restrictive than 



a smallness condition on the discretisation of the L ? norm of po (see Section]^. 

Remark 4.2 Whenever an implicit Euler scheme w.r.t. the time variable t is 
applied to the discrete dynamical system ( |4.4l), o ne follow down to the space 
discrete Jordan-Kinderlehrer-Otto scheme |ll Il8j. analyzed in [7] in the limit 
case 7 = 0. 
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5 Blow-up criteria for the particle system 



The main idea to derive blow-up criteria for the dynamical system ( |4.4[ ) is to 
follow carefully the time evolution of the square of the euclidean norm |Af(t)|, 
giving the approximation of the second moment of p by (3.6). Due to the 
gradient structure (4.4 1, it holds 

Moreover, since 

G y [XX} = —h(N — 1) log A + U[X] - A~ 7 W 7 [X] , 

we have finally 

^l|Af(t)| a = / l (JV-l)-7W r [*](«) = ft(JV-l)+7Gf 7 [^](t)-7Cr[Af](t). (5.1) 

It remains to find lower bounds for W 7 and U, to be plugged in the differential 
equation (5.1). To begin with, we denote hereafter A the (N — 1) x N matrix 
of the system Yi = X i+1 — Xi, i — 1, . . . , N — 1, B the N x TV matrix of the 
system obtained adding to the previous system the equation Yli=i Xi = 0, and 
we establish the following technical Lemma, giving the euclidean norm \X\ in 
term of the relative differences (Xj — Xi), j > i. 

Lemma 5.1. For any X E R N lying in the hyperplane X^iLi Xi = 0, it holds 

N-l N 

N\X\ 2 =Y1 E ( X j- X i) 2 - ( 5 - 2 ) 

i—1 j=i + l 

Moreover, if X E U := {X E R N : X, < X. l+1 and J2^Li Xi = 0} and 
Y = AX, it holds 

\X\ 2 =Y T (AA T )- 1 Y, (5.3) 

and 

N-l N-l 

W 2 ^^^^ 1 "* 1 "^' (5 - 4) 

»=1 j=i 

Proof. Observing that N Xf = Y^j=i(Xi — Xj)Xi, we have 

JV i-l N-l N 

N \x\ 2 = E E( x * - x i) x * + E E ^ - x ^ 

j=2 j=l i=l j—i+1 

N-l N N-l N 

= E E - x ^ + E E ( x * - x i) x i . 

i—1 j=i-\~X i—1 j—i-\-l 
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and the identity (^ follows. Next, with Y = (Y, 0) = (AX, 0) £ R N , we have 
\X\ 2 = IB- 1 ?] 2 = Y T (B B T )^ 1 Y = Y T (AA T )- 1 Y , 



and (5.3) is proved. Finally, observing that the matrix (AA T ) _1 is symmetric 
with positive entries each of which, outside the main diagonal is greater or equal 
to jj, and on the main diagonal greater or equal j?, we obtain inequality (5.4). 

□ 

Proposition 5.2 (Blow-up criterion induced by W 7 ). Let 7 £ (0, 1) and N > 3. 

Assume that the initial data X(0) £ il satisfies 



l*(0)| 2 < 



M\ ' +1 N 2 ^(N- 1) 
~2J (N + l)^ +1 



(5.5) 



Then, the corresponding solution X{t) £ IIo vanishes in finite time. Moreover, 
criterion (5.5) is incompatible with condition (4.6). 

Proof. From the finite form of the Jensen inequality applied to the convex func- 
tion (j)(z) = z~ 7 / 2 , we obtain for any X £ II 



(5.6) 



where /i(iV) 1 = N<,N 2 1 - > . Plugging (5.2) into (5.6), we obtain the lower bound 



W y [X] > 



1 



7 fx(N) \N-l 
and the differential inequality for \X(t)\ 2 



-7/2 



dt 



\X(t)\ 2 < 2(N - 1) — hN(N - 1) 



\X\ 



N - 1 



(5.7) 



7/2 



\x(t)\- 



Criterion (5.5) follows immediately. The incompatibility with condition (4.6), 
is a direct consequence of the right inequality in (4.12) and of the lower bound 
for the discrete interaction energy 



W^X] >-(N-l), 



obtained plugging (5.5) into (5.7). 



□ 



It is worth noticing that the lower bound (5.7) is the approximation of (2.7) 
and that the former has been obtained exactly as the latter, i.e. using the 
Jensen inequality applied to the same convex function. Moreover, the discrete 
criterion (5.5) converges toward the continuous criterion (2.9) as h — > 0, if the 
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initial data X(0) G Hq is chosen such that ft \X(0)\ 2 — > J H |aj| 2 po(20 dx as ft — » 0. 



It is much more tricky to derive lower bounds for the discrete entropy (4.1 1 
giving the corresponding blow-up criteria. The reason is that one has to use the 
finite form of the Jensen inequality in a sharp way, as it has been done in (5.6). 
Let us first deduce an intermediate criterion to illustrate this technical problem. 

Proposition 5.3 (Blow-up criterion induced by U). Let 7 £ (0, 1) and N > 3. 

Assume that the initial data X(Q) G ITo satisfies 



\X(0)\ 2 < h 2 (N - 1) e- 2/ ~< cxp (- 



h(N-l) 



G 7 [*](0) 



(5.8) 



Then, the corresponding solution X(t) G Ho vanishes in finite time. Moreover, 
criterion (5.8) is incompatible with condition (4.6). 



Proof. Given X £ Ilo, from (5.4) and the finite form of the Jensen inequality 
applied to the concave function log z we derive 

N-l N-l 



logixr > log 



log M (iV) £ (X i+l - Xi)(X j+ i - X } ) 

N-l N-l 



> io g (jv - 1) + ^N) E lo s((^+i - ^)(^+i - x i)) 

z=l j=i 
N-l 

= log(JV - 1) + Nfi(N) J2 log(Jf<+i - *0 ■ 



Hence, being 17 [J£] = /i(JV - 1) log ft - ftE^ 1 lo g(^+i - X i), & nolds 



U[X] > h(N - 1) log ft 



(N-l) 



log\X\ 



(N-l) 



\og(N - 1) . (5.9) 



Plugging the lower bound (5.9) in (5.1) and using the decreasing behaviour of 
G 7 [<Y](f), we get the differential inequality 



dt 



\X(t)\ 2 < 2 (N - 1) + 2 7 ft- 1 G 7 [A'](0) + j(N - 1) log \X(t)f 
-j(N- l)log(ft 2 (iV- 1)) . 



Criterion (5.8) follows immediately. 



To prove that (5.8) is not compatible with condition (4.6), it suffices to write 



(5.8) equivalently as 



h(N-l) 



\o g \x(o)\ 2 +u[x](o)- 



h(N-l) 



\og(h 2 (N-l)) < W^[X](0)— (N-l) 



2 2 ' 7 

and to observe that the quantity in the l.h.s. of the above inequality is non- 
negative by (5.9). The claim follows then by using the right inequality in (4.12). 

□ 



17 



It is immediately seen that criterion ( |5.8[ ) doesn't converge toward (2.10) as 
h — > 0. This is essentially because the lower bound (5.9) has been obtained 
using the non-sharp inequality (5.4) ((5.4) becomes an identity iff N = 3). As 
a consequence, (5.9) is not consistent with (3.9). However, (5.8) is contained in 
a continuum of blow-up criteria obtained mimicking the proof giving (3.9). 

Proposition 5.4 (A continuum of blow-up criteria induced by U). Let'y £ (0,1) 
and N > 3. Assume that the initial data X(0) 6 IIo satisfies 



\X(0)\ Z < h z (N - iy C(N) exp ^- 

where C(N) := sup C'(fi), the supremum being taken over all fi = (/zo, . 
with Hi > 0, i = 1, . . . , N — 1, arbitrarily fixed, fiQ = /in = 0, and 



C( M ) 



,n y 1 ( 2 N ~ x \ 

Y^i - pti-i) 2 ) exp - Y lo SMi ] 

\»=1 / \ i=l / 



(5.10) 
(5.11) 



Then, the corresponding solution X(t) G IIo vanishes in finite time. Moreover, 
criterion (5.10) is incompatible with condition (4.6). 

Proof. With [i arbitrarily given as above and for lell, we have 

N-1 N-1 

U[X] = h(N - 1) log h - h Y log((X i+1 - XJm) + h Y log Mi 



'N-1 



> h(N - I) log h - h(N-l) lag I Y jfri ( - X ^ 1 ~ X ^) + h Y, lo %f 1 



N-1 



. i=l 



i=l 



N 



N-1 



= h(N - 1) log(h(N - 1)) - h(N - 1) log Y MtM-i - Mi) )+hY lo S^ 

\i=l / i=l 

h(N-l) 



> h(N - 1) log(/i(iV - 1)) - H( - N n ^ log(|A| 2 



logC(M)- 

(5.12) 



It remains to prove that the positive constant C(/i) is upper bounded, so that 
C(N) is finite. Indeed, setting Jl = (/ii, . . . , /!tv__i) t , from the one hand we have 



N 



N-1 



^(M,-M,-i) 2 =JI t (AA t )JI> Ai Y 



i=l 



i=l 



where Ai — 4 sin (ow) is the smallest eigenvalue of the squared matrix A A 
and from the other hand 



N-1 



N-1 



1 



N-1 



exp 



N — 



2 

i ' 



(=1 
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by the Jensen inequality. Hence, C(p) < (Xi(N — We conclude plugging 

the lower bound for U[X] obtained optimizing (5.12) w.r.t. /x, into (5.1) and 
using again the decreasing behaviour of G 7 [X] (t) , so that 

d 



dt 



\X(t)\ 2 <2(N-1) + 2 1 h- 1 G 1 [X](0) + 7 (JV - 1) log(\X(t)\ 2 ) 
- 2-f(N - 1) lag(h(N - 1)) - j(N - 1) log C(N) . 



Criterion (5.10) follows immediately. The incompatibility of ( 5.10[ ) with condi 



tion (4.6) is proved exactly as in Proposition 5.3 



□ 



It is actually possible to prove that the new lower bound (5.12) for the 
entropy U[X] is consistent with (3.9) and that criterion (5.10) converges toward 

(5.13) 



(2.10) as N goes to infinity, since the constant C(N) in (5.10) satisfies 

1 



lim N- 1 C(N) 
n^+oo 2ne 



Corollary 5.5. The limit (5.13) holds true. 

Proof. First we assume M = 1, without loss of generality, and we prove that 
lim sup logO^- 1 C(N)) < -log(27re) . (5.14) 



The latter is an immediate consequence of (5.12). Indeed, noting H the inverse 
of the distribution function G of -J=e~ x / 2 , for X = (H ( N l + i ))i= i,. ..,iv £ n and 
appropriate 9i £ [0, 1], we have 
2 



log(AT- 1 C(N)) < 



h(N-l) 



U[X) - 2 \og(h(N - 1)) + \og{\XY) - logiV 



N-l 



(N 
log 



T)£ log r 



N+ 1 



1 



N 



N - 

l0g(27T) 



N+l 

JV-l 



-21og(MJV-l)) 
N+l' 



log 



log 



1 



N - 



N 

T 5> 



N + l 



N+l 



log 



N 

- 2log(h(N - 1)) 

N + l' 
N 



Passing to the limit as N goes to infinity, we obtain (5.14). The limits above 
can be easily justified using the increasing behavior of H and its symmetry with 

respect to m = |, so that J Q H 2 (m)dm = 2 J Q 2 H 2 {m)dm = 2 jl H 2 (m)dm. 

For instance, if N = 2n + 1, we have 



I 



1 



N 



N 



N - 



N+l 



N ■ 



N + l 



19 



H 2 (m) dm < I N < 2 / H 2 (m)dm 



and 



To conclude, it remains to prove that there exists /i = (0, fi\, . . . , [1n-i, 0), 
with (if > 0, i = 1, . . . , N — 1, such that TV -1 C(fi) converges toward (27re) _1 
as N goes to infinity. So, we consider /i, = G'(H(jj)), i = 1, . . . , N — 1. Then, 
from the one hand we have, for i = 1, . . . ,N — 2 and appropriate 6i £ [0, 1], 



so that 

N 



W-2 



H' 



i + e 

N 



1 H l + ° 



N 

i=l i=l 

From the other hand 

2 

/N-l \ N-l 



N 



(5.15) 



JV-l 



n 

i=l 



Z7T 



TV 



Therefore, using (5.15) and (5.16), we obtain the limit and the proof. 



(5.16) 



□ 



Finally, we shall show that the non optimal family of criteria obtained from 



(5.10) when C(TV) is replaced by C(/i), contains criterion ( |5.8[ ). The latter is 
verified obviously iff C(fi) = ]y— j. With the constant C(/i) in (5.11) written 
equivalently as 



N-l 



c(m) = 2 ~ fo-i) 



'N-l \ w^T 



and /ii = > 0, /^i = 1, i = 2, . . . , TV — 1, the identity C(/i) = becomes 



W- 1 = (v - v + 1) , 

N-V ' 



and the above equation has at least one solution for any TV > 3, (exactly one 
solution for TV = 3 and two for TV > 3). 



Remark 5.6 It is easily seen that the blow-up criteria (5.5), (5.8) and (5.10) 



are all invariant w.r.t. the scaling (h\, X\) — (A 7 h, XX) of the discrete setting. 



6 The discrete threshold phenomenon for 7 = 

In the limit case 7 = 0, the interaction potential kernel if 7 has to be replaced 



with Kq(x) — — log |x| in equation (1.1). This problem reproduces in any di- 
mension d > 1, the well know critical mass phenomenon of the two dimensional 
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parabolic-elliptic Keller-Segel system [5] , basically because both the internal en- 
ergy U(p) and the interaction energy Wo(p) scale additively with a logarithmic 
correction, under the action of the mass-preserving dilation f x (x) = jg /(f)) 
i.e. 

£a[p X ]=M fy-d) log\ + £ [p}. 
Then, the energy is bounded from below on the set of integrable densities p with 



finite energy iff M < M c = 2d, (see also (2.5)) 



This problem has been analyzed, for any dimension d > 1, in the Eulerian 



formulation ( 1.1 ) in [TTj . and following the gradient flow formulation (3.8 1, writ- 
ten in self-similar variables, in [7J (see also [5]). More specifically, the authors in 
[7] proved the convergence of the Jordan-Kinderlehrer-Otto scheme [IS] adapted 
to that problem when M < M c . However, without going through this generality, 
if d = 1 it is still possible to define the particle scheme when 7 = 0, setting 

N-l N 

W [X] := -h 2 log(Xj - Xi) and X(t) = -VG [X](t) . (6.1) 

i— 1 j=i-\-l 

Then, thanks to the homogeneity property of the corresponding discrete energy, 
i.e. 

h 2 

G [AX] = -h(N - 1) log A + — N(N - 1) log A + G [X] , 

it holds 

^WI'-^-i)(i-^ 

and the solution blows up iff h > jt, i.e. iff 

M > 2 • (6-2) 

It is worth noticing that the above discrete blow-up criterion converges as 
N — !• oo toward the continuous blow-up criterion M > 2. 

We conclude this short review of the case 7 = with the proof of the global 
existence of a solution of the particle system under the condition reverse to 



(6.2), showing that the threshold phenomenon holds true also at the discrete 



level. To the best of our knowledge, this is a new result. 
Theorem 6.1 (Global existence). Let 7 = and N > 3. If 

M < 2 , (6-3) 

then for any <-f(0) € II there exists a unique global solution X{t) G II of the 
dynamical system (6.1). Moreover, if N = 3, the discrete entropy U[X\(t) is 
decreasing in time. 
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Proof. The proof of the local existence of a trajectory X(t) starting from X(0) 
is exactly the same as for 7 £ (0, 1). To prove that this trajectory is global in 
time under condition (6.3 1, i.e. none of the contiguous particles collide, it is 



enough to prove that the discrete entropy U[X] stay upper bounded along the 

trajectory X(t) globally in time, since U[X] = ft log (jjf^ 1 ft pQ+i — X;) -1 ^ 

goes to +00 as a relative difference (Xi+\ — Xi) vanishes. 

The proof is inspired by the theory of the two dimensional Keller-Segel sys- 
tem. Indeed it is based on the following (non-optimal) discretisation of the well 
known logarithmic Hardy-Littlewood-Sobolev inequality |12j . 



2W [X] < hNU[X] - h 2 N(N - 1) log ft, X £ U, 



(6.4) 



that can be obtained exactly as the right inequality in (4.12). Let us observe 



that the non-optimality concerns only the constant in the r.h.s. of (6.4) that 
goes to +00 as ft goes to 0, and not the constant factor of J7[X], That is 
the reason why we are nonetheless able to obtain the sharp global existence 



condition (6.3). So, as in the continuous setting, using (6.4) and the decreasing 



behavior of the discrete energy Go[X] along the trajectories, we have for 9 > 

(1 - 9)U[X](t) < G o [X](0) - eu[x](t) + W [X](t) 

2d 



<G [X}(0)+[ 1 



hN 



W [X](t) -0h{N- 1) log ft. 



Choosing 6 



hN 



, we obtain a uniform in time upper bound for f/[Vf](t) under 



condition (6.3). 



Finally, if N = 3, the decreasing behavior of U[Af](i) relies on the decreasing 



behavior of the function <j>o(X) := —ft X^ili 1 l°g(^i+i — limit of the function 
4>-y(X) as 7 goes to 0. Indeed, <fio(X(t)) evolves according to the equation 

A M*(t)) = 2ft 2 (ft - 1) {{x 2 - x x y 2 + {x 3 - x 2 )- 2 - (x 2 - x 1 )- 1 (x 3 - x^- 1 ) 



dt 



+ h 3 {x 2 -x 1 )-\x 3 -x 2 )- 1 



Since ft < 3 by (6.3 ), applying the Jensen inequality the above equation becomes 



dt 



M*(t)) < ft 2 (3ft - 2){X 2 - Xi)" 1 ^ - X 2 )~ l 



and the claim follows. 



□ 
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Figure 1: Behaviour of the three particles scheme, with M = 1 and 7 = |, in the 
phase plane (u, w) = ((X2 — -Xi)> (X3 — ^2))- Top Ze/t : the blow-up criterion (5.5 1 
is figured as a bold line. Top right : the blow-up criterion ( 5. 10 1 is figured as a bold 
line. Middle left: the curve 7 W 7 [X] = h(N — 1) of the critical points of the energy, 
is plotted in bold. Middle right : the global existence criterion (4.61 is figured as a 
bold line. Bottom left : the unstable manifold starting from the energy maximal point 
and separating the two basins of attraction (global existence and blow-up) has been 
computed numerically and plotted in bold. Bottom right : all the previous lines are 
plotted on the same figure for the sake of comparison. 
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